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Abstract 

We study dephasing in semiconductor-superconductor structures caused 
by coupling to a voltage probe. We consider structures where the semicon- 
ductor consists of two scattering regions between which partial dephasing 
is possible. As a particular example we consider a situation with a double- 
barrier junction in the normal region. For a single-mode system we study 
the conductance both as a function of the position of the Fermi level 
and as a function of the barrier transparency. At resonance, where the 
double-barrier is fully transparent, we study the suppression of the ideal 
factor-of-two enhancement of the conductance when a finite coupling to 
the voltage probe is taken into account. 
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1. Introduction 

Charge transport through a normal conductor-superconductor (NS) interface is 
accompanied by a conversion of quasiparticle current to a supercurrent. In the 
Andreev reflection, by which the conversion occurs, an electron-like quasiparticle 
in the normal conductor (with an excitation energy lower than the energy gap 
of the superconductor) incident on the NS interface is retroreflected into a hole- 
like quasiparticle (with reversal of its momentum and its energy relative to the 
Fermi level) and a Cooper pair is added to the condensate of the superconductor 
0. For an ideal NS interface, a signature of Cooper pair transport and the 
Andreev scattering is a doubling of the conductance compared to the normal 
state conductance. 
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A theoretical framework for studying scattering at NS interfaces is provided 
by the Bogoliubov-de Gennes (BdG) formalism where the scattering states 
are eigenf unctions of the BdG equation which is a Schrodinger-like equation in 
electron-hole space (Nambu space). 

For a phase-coherent normal region, the conducting properties can be found 
using a scattering matrix approach based on the BdG equation. In the linear- 
response regime in zero magnetic field, Beenakker || found that the sub-gap 
conductance G = dl/dV is given by 

2 N rp 2 

G m = 2G Tc (tttp-ttt]- \ =2G V " (1) 

n=l \ Z 1 n) 

which, in contrast to the Landauer formula ||, [|, ||, [7j for the normal state 
conductance, 

TV 

G N = G Trttt = Go ^T n , (2) 

71=1 

is a non-linear function of the transmission eigenvalues T n (n = 1, 2, . . . , N) of 
tv . Here Go — 2e 2 //i is the quantum unit of conductance for a spin-degenerate 
system and t is the N x N transmission matrix of the normal region, N being 
the number of propagating modes. The conductance formula holds for an arbi- 
trary disorder potential and is a multi-channel generalization of a conductance 
formula first obtained by Blonder, Tinkham and Klapwijk || who considered a 
delta function potential as a model for the interface barrier potential. Eq. (|]) 
is computationally convenient since standard methods developed for quantum 
transport in normal conducting semiconductor structures can be applied to cal- 
culate the transmission matrix t or the corresponding transmission eigenvalues 
T 

Equations (0) and (||) apply to two-probe structures where the length L of the 
disordered normal region is much shorter than the phase-correlation length over 
which quasiparticle propagation remains phase-coherent. 

The effect of coupling a mesoscopic region to a voltage probe was first used by 
Btittiker in the context of the role of quantum phase-coherence in normal con- 
ducting series resistors and subsequently in the study of the cross-over from 
coherent to sequential tunneling in series of barriers [[H]]. In the work of Chang 



and Bagwell [|Tl|] the model of Biittiker [JTOj was used as a normal-metal probe 
controlling the Andreev-level occupation in a Josephson junction. In this paper 
we apply the model of Biittiker to a two-probe semiconductor-superconductor 
junction with a phase-coherent normal region. However, here we consider the 
case where the dephasing reservoir is constituted by a real voltage probe. As a 
general model we consider a normal region with two arbitrary scattering regions 
(characterized by scattering matrices Si and S 2 ) separated by a coherent bal- 
listic conductor of length L to which the voltage probe couples, see (a) in Fig. 

B 
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Figure 1: Scattering scheme for a two probe structure with two scattering regions (charac- 
terized by the S-matrices Si and S2 ) between which dephasing is caused by coupling to a 
voltage probe. The left probe (1) is in the normal state and the right probe (2) is a metal 
which can be either in the normal or superconducting state. The voltage probe could e.g. be 
an STM tip. Panel (a) shows the generic structure, (b) corresponds to the present model for 
the voltage probe, (c) to the phase-coherent regime with no coupling to the voltage probe, 
and (d) to the sequential tunneling regime with full coupling to the voltage probe. The arrows 
indicate possible directions of current flow and the arrows with dotted midsections indicate 
equilibration of quasiparticles. 



Specifically we study the case where the total current is conserved, but, due to 
the coupling to the voltage probe the non-equilibrium distribution of electrons 
and holes equilibrates to a Fermi-Dirac distribution. In Nambu space this cor- 
responds to a decay of electron-hole states in the normal region. As a particular 
example we consider a semiconductor-superconductor structure with a normal 
region containing a double-barrier junction of length L and model the cross-over 
from vanishing dephasing to full dephasing in the conducting properties. 

We note that in our approach the voltage probe is real as it is also the case 
in the recent work of Gramespacher and Buttiker where a tunneling tip is 
coupled weakly to the normal side of a two-probe NS structure. As in the original 
work of Buttiker |9|, |10j the voltage probe could also be considered as a fictitious 



probe used in modeling a finite dephasing length of the normal region. However, 
we note that the model only accounts for dephasing caused by inelastic scatter- 
ing and that thermal dephasing, which is the dominant dephasing mechanism 
in most semiconductor-superconductor structures, can not be accounted for by 
this model |13|]. Dephasing due to a finite temperature can in principal be in- 



corporated by calculating the full energy dependence of the S-matrices, see e.g. 



14] , |15l , jig , [Tfll , or through Green functions methods, see e.g. |17]] and references 
therein. 

The paper is organized as follows: In Section II the S-matrix formalism is 
introduced, in Section III we formulate our general model and present the cor- 
responding scattering scheme. In Section IV we present results of an application 
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of our scattering scheme to the problem of double-barrier tunneling. Finally, in 
Section V discussion and conclusions are given. 



2. Scattering matrix formalism 

Using the Biittiker voltage probe [[HJ we need to consider four-probe structures 
(two real and two probes accounting for the voltage probe) and we therefore 
need multi-probe generalizations of Eqs. ([[]) and (g). The scattering approach to 
dc transport in superconducting hybrids follows closely the scattering theory de- 
veloped for non-superconducting mesoscopic structures, see e.g. the monograph 
of Datta Jl8| . For a multi-probe system, the normal state current is given by the 
multi-probe conductance formula of Biittiker pj which in the linear-response 
limit gives the following current in lead p: 



M 

X V = Wp ~ V q\ > Gpq = G Tl 5 P A' ( 3 ) 

where q = 1, 2, 3, . . . , M labels the M probes and V q is the potential of probe 
q. The S-matrix S is an M x M block-matrix with matrices S pq describing 
scattering of an electron state in lead q into an electron state in lead p. These 
matrices are N p x N q matrices where N p and N q are the number of modes in lead 
p and q, respectively. In the linear-response limit, the S-matrices are evaluated 
at the Fermi level. 

For an NS interface with multiple normal probes, Lambert, Hui and Robinson 
fli~4| considered the sub-gap current which in the linear-response regime is given 
by the formula 



M 

Ip — ^ G pq [V g - V] , G pq = G Ti (l5 pq + S p h q {S C p q Y - S c pq {S pq } f J , (4) 

q=l 

where 1 is an N p x N p identity matrix, S pq is the Kronecker delta symbol and 
V is the potential of the superconductor. The S-matrix in electron-hole space is 
given by 



(gee gch \ 
S hc S hh ) ' (5) 

where the sub-matrix S pq {S pq ) is for normal scattering of an electron (hole) 
state in lead q to an electron (hole) state in lead p, and S pq (St%) is for Andreev 
scattering of an electron (hole) state in lead q to a hole (electron) state in lead 
p. 

The two-probe case (M = 2) was studied by Takane and Ebisawa [19] and 
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subsequently by Beenakker || who derived Eq. ([]]) within the Andreev approx- 
imation 0] and the rigid boundary condition for the pairing potential of the 
superconductor. Within these approximations, the interface acts as a phase- 
conjugating mirror and the interface scattering can be described by an S-matrix 
which at the Fermi level is given by [TJ, [| 



I e^Jl \ 

where tp is the phase of the pair potential of the superconductor. 

The form of the BdG equation || in the normal region where the pairing 
potential A is zero means that the elastic scattering in the normal region due to 
an arbitrary disorder potential is characterized by the block-diagonal S-matrix 



where S is the normal state S-matrix for electrons (the S-matrix entering e.g. Eq. 
(§) ) and S* is the corresponding S-matrix for scattering of holes. The S-matrix 
in Eq. (|) entering Eq. (f|) is given by the composite result 5* = Sa <8> <Sn where 
we use the notation of Datta fl8| with the meaning of ® found by elimination 
of internal current amplitudes. 

3. Scattering scheme 

We consider an NS structure with a normal region containing two scattering 
regions 



connected by a ballistic conductor of length L. The free propagation of electron- 
like quasiparticles is described by 

with elements given by X nn i = 5 nn i exp (ik n L) . For a wire of width W with a 
hard-wall confining potential the longitudinal wave vector k n at the Fermi level 

is given by k n = k-p\jl — (rm/k-pW) 2 . 

We connect voltage probe to the center of the conductor, see panel (b) in Fig. 
|T]. Using the model of Biittiker [US] the conductor is now characterized by the 
unitary matrix 
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(10) 



/ 



which is a generalization of the S-matrix in Ref. 0] with the matrix X taking 
the finite length of the conductor into account. Here < ( < 1 is a parameter 
determining the coupling of the conductor to the voltage probe . The scattering 
from lead q to lead p is given by the sub- matrices {S^} which are zero along 
the diagonal since there is no back-scattering from the conductor and the other 
zeros reflect that lead 1 couples only to leads 2 and 3, lead 2 only to leads 1 
and 4 etc. The basic idea is that current can flow either directly through the 
phase-coherent conductor, as shown in panel (c) of Fig. [I], or in the sequential 
way, as in panel (d). Using the scheme in panel (b), the cross-over from coherent 
to sequential tunneling can be modeled by changing the coupling from ( = 
to ( = 1. In this approach, probes 3 and 4 have a common reservoir (that of 
the voltage probe) and quasiparticles entering probe 3 are re-injected through 
probe 4, and vice versa, so that the common reservoir of probes 3 and 4 does 
not supply or draw any net current (7 3 + J 4 = 0). 

The composite S-matrix of the normal region S = S 2 <S> <g> S\ can now be 
found. Eliminating the internal current amplitudes and introducing the matrices 
7i 2 = [i - (1 - C) Xr 2 Xr[] _1 and 721 = [i - (1 - C) Xr[Xr 2 ] we get 



S 



( n + (1 - C) t' 1 ^i 2 Xr 2 Xt 1 VT^Ct' lll2 xt' 2 

VT=(t 2l21 Xti r' 2 + (1 - C) t 2l2l Xr\Xt' 2 

vcx i/2 ^ 2 h vca-c)* 172 ^^ 

V v /C(l-C)^ 1/2 r 2 72iXt 1 VCX 1 / 2 ^ 



VCAl nX 1 / 2 VC (1 - Ot'iluXr.X 1 / 2 \ 

VCii-C^iXrix 1 / 2 VCh^x 1 / 2 

CX^r'^X 1 / 2 -;/R{i - (X^r'^nXr.X 1 / 2 } 

-VW {1 - C^ 1/2 r 2 72iXr' 1 X 1 /2} cx l/2 r272lX i/2 J 

(li; 

which together with Eq. (§) gives the normal state conducting properties. 

In order to apply Eq. (^) for the conducting properties of the superconducting 
state we need to calculate the two block-matrices S cc and S hc of the composite 
S-matrix S = Sa <S> Sn- Since the phase (p of the pairing potential for a 2DEG-S 
system (in contrast to S-2DEG-S systems) can be arbitrary and the conducting 
properties are independent of (p, we may for simplicity let (p = 0. Eliminating 
internal current amplitudes we get 
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5i 3 

6 

o s 33 



iSi252 2 r 22 iS23 

6 

£32<!>22r22'S23 



s h 



( S\\ — <Si2<?2 2 r22 l -'2l S^^^ — S-\iS%£ s 

6 

^31 — >^32'522r22'521 
\ t>41 — <S42i?22r22<521 ^43 — ^A2^22^ 22^23 
/ —i<Sl2^22<S21 — iiSi 2 r22<S23 — 2iSi 2 r 22 iS 2 4 

6 6 6 6 

— iiS3* 2 r 22 iS2l — ^iSjj 2 r22<$23 _ 2iS| 2 F 22 iS 2 4 
\ — iS^2^22^21 — iSl2^22^23 ~i^A 2 ^22^2A / 



11 

^34 
1S44 

\ 



iSi25 22 r 2 2iS24 \ 

6 

l ^32'522r22'524 

£42 i5| 2 r 22 524 / 



2.) 



(13) 



where the matrix 1^22 = [l + <S 22 <S 22 \ has been introduced. The zero matrices 
in row 2 and column 2 reflect that the superconductor (probe 2) does not carry 
any quasiparticle current due to the gap in the density of states. At the interface, 
the quasiparticle current is transformed to a supercurrent carried by the Cooper 
pairs of the condensate so that current is conserved. For a finite bias exceeding 
the energy gap of the superconductor, the current in the superconducting probe 
will be carried by both quasiparticles and Cooper pairs. 

To find the normal state conductance Gn we apply Eq. (|3]). From Kirchhoff's 
law Yl 3 P = and the condition Z 3 + Z 4 = for the voltage probe get 



Xi 




Gu + Gn + Gu —G13 - Gu 

-031 - Gil 031 + 032 + 041 + 042 



Vi-V 2 

v-v 2 



(14) 



where V — V3 = V4 is the common potential of the voltage probe reservoir. The 
conductance Gn = dX 1 /d{Vi — V 2 ) is therefore given by 



Gn = G12 + G 



(013 + 014) (031 + G 



41 ) 



13 



'14 



'31 



032 + 041 + 



12 



(15) 



In the superconducting state we apply Eq. 
that 



and in the same way we find 



h 




Gu 
G31 + G41 



G13 + G14 

G33 + G34 + G43 + G44 



V 1 -V 2 

V -v 2 



(16) 



so that the conductance G 



NS 



G 



NS 



Gu — 



dIi/d(Vi — V 2 ) is given by 

(G13 + G14) (G31 + G41) 



G33 + G34 + G43 + G 



(17) 



t 11 



Despite the minus signs in Eqs. (|15|) and (|i"7[), both Gn and Gns are positive 
definite which follows from the unitarity, time- reversal, and electron hole sym- 
metries of the S-matrices f2Cj. Eqs. ( |15"1) and (|17D together with Eqs. (|Tl"D , (|i~2"D, 
and ( |T3"D form the basis for our calculations of the conducting properties of two- 
probe semiconductor-superconductor junctions with a normal region containing 
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two scattering regions separated by a phase- coherent region which couples to a 
voltage probe. 

The sequential tunneling limit (£ = 1) corresponds to series connected resis- 
tors, see panel (d) in Fig. [I]. Thus we can in general combine Eqs. ([!]) and (||) to 
obtain 

Gn ^= G ( h „ — ) , Gns ^— G 



Erp 1 \ ^ rp I 5 IN o U I \— ^ rp 1 



T 2 



(2-T 2 ,n) 2 



where T n l and T n 2 are the transmission eigenvalues of t\t\ and t 2 4' respectively. 

The conductance in the zero-coupling limit (( = 0) can of course be calculated 
directly from Eqs. ([[J) and (g) with the transmission eigenvalues T n i 2 of £12^12 
where t\2 is the transmission matrix of the composite S-matrix S12 = S2®U®Si 
see e.g. the structure of Eq. (181). 



12- 
■ 



4. Double-barrier tunneling 

As a specific example we consider double-barrier tunneling with two identical 
barriers which we model by delta function potentials V(r) = H [5(z) + 5(z + L)]. 
Here L is the separation of the two barriers and H is the strength of the poten- 
tials. The S-matrices are given by 

c c ( r s t s \ ( , 1 , , -iZ/cos6 n 

S 1 =S 2 = . „ , (*,)__, = 5 nn , , [r s ). = 6,, 



t s r s J ' v " nn 1 + iZ/cos9 n ' v Jnn 1 + iZ/ cosOn 

(19) 

where Z = H/hvp is a normalized barrier strength and cos# n = k n /kp takes the 
transverse momentum of mode n into account [EH], E31. 



In the phase-coherent limit we calculate the transmission matrix t 12 corre- 
sponding to the composite S-matrix <Si2 = S2 <8> C/ <8> «Si. Since the transmission 
is diagonal, £12^12 nas transmission eigenvalues given by 



T 12 , n = (l + 2Z e 2 (0 n ) [1 + costf n + 2Z e (9 n ) sintf n + Z c 2 (fl n ) (1 - costf n ) ]) , 

(20) 

where Z e (6 n ) = Z/ cos8 n and = 2kpL cos8 n . Comparing to the transmission 
eigenvalues Ti )n = T 2jn = [1 + Z%(9 n )] 1 of the individual barriers it is seen 
that for the double-barrier system, T 12) n can be equal to unity even for a finite 
value of Z which is in contrast to the result for Ti jra = T 2 ,n- This is one of 
the important differences between phase-coherent quantum transmission and 
classical transmission probabilities. At the values of Z where mode n becomes 
fully transparent a peak in the conductance is to be expected. Apart from the 




Figure 2: Double-barrier junction with a normalized barrier strength Z = 0.05 and an aspect 
ratio L/W = 3. Panel (c) shows the normal state conductance Gn, (b) the conductance Gns 
of the superconducting state, and (a) the normalized conductance Gns/Gn for different values 
of the dephasing parameter £. The plots are shown as a function of k?W '/it < 2 corresponding 
to only a single propagating mode. Notice the different scales on the conductance axes. 
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trivial solution Z = we find that Z = — cos8 n cot($ n /2) which means that the 
mode can be fully transparent for 7r < -d n < 2tt (modulo 2n) only. 

As a numerical example we consider a double-barrier junction with an as- 
pect ratio L/W = 3. Fig. || shows the conductance as a function of kpW/7r 
for Z = 0.05. For k-pW/ir < 1 there are no propagating modes and for 1 < 
k^W/n < 2 a single mode is propagating. In panel (c), the normal state con- 
ductance Gn increases from zero to approximately Go at the onset of the first 
propagating mode. As the coupling parameter is increased from towards 1, 
Gn approaches G /2 and it actually becomes a little lower than G /2 due to 
the finite barrier strength (at k F W/n = 2, Eq. ([18]) gives G N /G = 150/301). 
In the coherent regime, oscillations caused by size-resonances are seen and as 
the sequential tunneling regime is approached these resonances vanish. In panel 
(b), the conductance Gns m the superconducting state shows the same over- 
all behavior as seen for Gn, even though the size-resonances are much more 
pronounced. These size-resonances are even more pronounced in the normalized 
conductance Gns/Gn shown in panel (a). As the coupling parameter is increased 
the resonances vanish and the normalized conductance is suppressed below the 
ideal factor-of-two for the enhancement of Gns compared to Gn- For these pa- 
rameters, ideal transmission and maxima in the conductance are expected for 
k F W/ir ~ 1.018,1.113,1.296,1.531,1.808 which is in full agreement with the 
plots in Fig. |2] where these numbers are indicated by arrows. 

In Fig. H we consider the Z-dependence of the conductance for a Fermi wave 
vector corresponding to kpW/ir = 1.9 so that only one mode is propagating. 
Panels (b) and (c) show the conductance Gns in the superconducting state and 
the normal state conductance Gn, respectively, as a function of the normalized 
barrier strength Z for different values of the coupling parameter (. For these pa- 
rameters, a peak in the conductance is expected for Z ~ 1.626 which is indeed 
what is clearly seen in the coherent tunneling limit (£ = 0) where this number 
is indicated by an arrow on the Z-axis. As the sequential tunneling regime is 
approached (£ — > 1) the peak is suppressed and an overall lowering of the conduc- 
tance is seen. However, dephasing may also suppress a destructive interference 
in the coherent transmission resulting in an increased conductance. In this par- 
ticular case this is seen for Z > 2.2 in both Gn and Gns and also in a range 
around Z ~ 0.7 in Gns where the conductance curves are crossing. Panel (a) 
shows the corresponding normalized conductance Gns/Gn which has the value 
of two for ideal Andreev scattering and decreases below the ideal factor-of-two 
enhancement when Z is increased so that normal scattering becomes dominant. 

For C = 1 and Z = Eq. (0) gives that G N /G = 1/2, G NS /G = 2/3, and 
Gns/Gn = 4/3. These numbers are born out by the numerical calculations as 
indicated by arrows on the conductance axes in Fig. ^ and they can also be seen 
in Fig. Q where there is a very low but finite barrier strength. 

In experiments, the barrier strength will often not be a tunable parameter. 
However, for a given barrier strength the conductance can be lowered or en- 
hanced due to phase- coherent quantum interferences in the double-barrier junc- 




Figure 3: Double-barrier junction with an aspect ratio L/W = 3 and a Fermi wave vector 
corresponding to kpW/n — 1.9. Panel (c) shows the normal state conductance Gn, (b) the 
conductance Gns m the superconducting state, and (c) the normalized conductance Gns/Gn 
as a function of the normalized barrier strength Z . The dephasing parameter £ is varied from 
to 1 in steps of 0.05. Notice the different scales on the conductance axes. 
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tion. Thus, increasing the coupling to the voltage probe, the dephasing can give 
rise to both an increasing or decreasing conductance depending on the specific 
parameters of the double-barrier junction. 

5. Discussion and conclusion 

In semiconductor structures and also semiconductor-superconductor junctions 
the conducting properties depend strongly on whether the length scale L of the 
device is much small than the phase-correlation length L c or it is comparable 
to L c . Another dephasing mechanism could be a coupling to a voltage probe. In 
this paper we have studied the effect of coupling a voltage probe to two-probe 
semiconductor-superconductor structures with two scattering potentials in the 
normal region. 

As a particular example we have studied the cross-over from phase-coherent to 
sequential tunneling in a double-barrier semiconductor-superconductor junction. 
For a finite coupling to the voltage probe, we find that size-resonances of course 
become less pronounced. This may increase or decrease the conductance when 
increasing the temperature depending on the specific parameters of the double- 
barrier junction. However, we also find an overall lowering of the conductance 
and a suppression of the ideal factor-of-two enhancement of the conductance 
compared to the normal state conductance. 
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